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Abstract 



We present limits on anomalous WW Z and WW^ couplings from a search 
for WW and WZ production in pp collisions at y/s = 1.8 TeV. We use 
pp — > evjjX events recorded with the D0 detector at the Fermilab Tevatron 
Collider during the 1992-1995 run. The data sample corresponds to an inte- 
grated luminosity of 96.0 ±5.1 pb _1 . Assuming identical WW Z and WWj 
coupling parameters, the 95% CL limits on the CP-conserving couplings are 
-0.33 < A < 0.36 (Ak = 0) and -0.43 < Ak< 0.59 (A = 0), for a form factor 
scale A = 2.0 TeV. Limits based on other assumptions are also presented. 

Submitted to Phys. Rev. Lett. 
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The vector boson trilinear couplings predicted by the non-Abelian gauge symmetry of 
the Standard Model (SM) can be measured directly in pair production processes such as 
qq — > W + W~, W ±r y, Zj, and W ± Z. Deviations from the SM couplings would signal 
new physics. Studies of such effects have been reported by the UA2 |], CDF |fj] and 
D0 collaborations. In this letter we report on the measurement of WWV couplings 
(where V — 7 or Z) using the diboson production processes pp —>■ WWX — > eujjX and 
pp — > WZX — > evjjX, where j represents a jet. 

The Lorentz invariant Lagrangian which describes the WW'j and WW Z interactions has 
fourteen independent coupling parameters |J, seven describing the WWj vertex and seven 
for the WW Z vertex. Assuming electromagnetic gauge invariance and CP conservation, 
the number of parameters is reduced to five: gf , Kz, k 7 , Xz and A 7 . In the SM at tree level, 
the coupling parameters have the values Agf (= gf — 1) = 0, Akz(= k<z — 1) = A/t 7 (= 
k 7 — 1) = 0, \z = A 7 = 0. The SM cross sections for pp — > W + W~X and pp — > W ± ZX 
production at the Tevatron, at y/s = 1.8 TeV, are 9.5 pb and 2.7 pb J7| respectively. 

Non-SM values of the coupling parameters would result in an increase of the production 
cross section, especially for large values of the transverse momentum of the W boson (pt). 
Since tree level unitarity restricts the WWV couplings to their SM values at asymptotically 
high energies, each of the couplings must be modified by a form factor e.g. Xz(s) = Az/(1 + 
s/A 2 ) 2 , where s is the square of the invariant mass of the WW or WZ system and A is the 
form-factor scale. We have used A = 1.0, 1.5 and 2.0 TeV. 

The analysis reported here uses pp — »■ eujjX events recorded with the D0 detector 
during the 1992-1993 and 1993-1995 Fermilab Tevatron Collider runs at yfs = 1.8 TeV, 
corresponding to a total integrated luminosity of 96.0 ±5.1 pb -1 . The D0 detector and 
data collection system are described elsewhere 0. The basic elements of the trigger and 
reconstruction algorithms for jets, electrons, and neutrinos are given in Ref. j5|. The analysis 
of evjjX events from the 1992-1993 Tevatron Collider run (13.7 ± 0.7 pb -1 ) was reported 
previously ||. This letter focuses on the analysis of the 1993-1995 data set of 82.3±4.4 pb -1 
and gives the combined results for both analyses. Further details are available in Ref. 0. 

The data sample was obtained with a trigger which required an isolated electromagnetic 
(EM) calorimeter cluster with transverse energy Et > 20 GeV and missing transverse en- 
ergy Ifjrp > 15 GeV. The offline event selection required that the EM cluster have \r]\ < 1.1 
in the central calorimeter or 1.5 < |?7| < 2.5 in an end calorimeter, where 77 is the pseudo- 
rapidity. Electrons were identified by requiring that the EM cluster pass the shower profile 
and tracking information criteria, as described in our earlier analysis ||. The presence 
of a neutrino was inferred from the ]ft T , calculated from the vector sum of the Ex mea- 
sured in each calorimeter tower. Jets were reconstructed using a cone algorithm with radius 



1Z = J (A77) 2 + (A0) 2 = 0.5. To remove spurious jets due to detector effects, this analy- 



sis used the same quality cuts as were used in Ref. [10|. Jets were required to be within 



\q\ < 2.5. The jet energies were corrected for effects of jet energy scale calibration, out-of- 



cone showering, energy from the underlying event [II]], and energy loss due to out-of-cone 
gluon radiation. 

The WW/WZ candidates were selected by searching for events containing an isolated 
electron with high E T (> 25 GeV), large ^ T (> 25 GeV) and at least two high E T jets 
(> 20 GeV). The transverse mass of the electron and neutrino system was required to be 
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consistent with a W -> ev decay (Mf = [2£^$ T (1 - cos0 eiy )] 1/2 > 40 GeV/c 2 , where <fr eu is 
the azimuthal angle between the electron and $ T vector). The invariant mass {m") of the 
two jet system (the largest invariant mass if there were more than two jets with E 3 T > 20 
GeV in the event) was required to be 50 < m" < 110 GeV/c 2 , as expected for a W — > jj 
or Z — > jj decay Monte Carlo studies showed that the dijet invariant mass resolution for 
signal events is 16 GeV/c 2 . The transverse momentum of the two gauge bosons was required 
to be within \pHp — p^\ < 40 GeV/c, as expected for WW/WZ production. 

There are two major sources of background to WW/ WZ — > evjj production: (i) W+ > 
2 jets with W — > ev; and (ii) QCD multijet events where one of the jets is misidentified as an 
electron and there is significant $ T in the event due to mismeasurement. Other backgrounds 
such as: tt production with subsequent decay to W + bW~b followed by W — > ev; WW or 
WZ production with W — > tv followed by r — > evv; ZX — > eeX, where one electron is 
mismeasured or not identified; and ZX — > ttX with r — > evv, are negligible. 

The W+ > 2 jets background was estimated using the vecbos |12[ event generator, 
with Q 2 = {pr) 2 -, followed by parton fragmentation using the herwig [TB| package and a 
detailed geant [14] based simulation of the detector. Normalization of the W+ > 2 jets 
background was determined by comparing the number of events expected from the vecbos 
estimate to the number of candidate events outside the dijet mass window, after the multijet 
backgrounds had been subtracted. The systematic uncertainties in this background are due 
to the normalization and to the jet energy scale correction. The multijet background was 
estimated following the same procedure used in our previous analysis ||. The background 
sample, which consisted of data events containing a jet satisfying the electron trigger se- 
lection but failing the electron identification, was normalized to the signal sample in the 
region $ T < 15 GeV where the actual WW/WZ contribution is negligible. The number 
of background events was then determined from this scaled background sample with the 
rest of the selection criteria applied 0. The backgrounds from tt — > W + bW~b and other 
minor sources were estimated using the isajet event generator [15] followed by detector 
simulation. Table I summarizes the background estimates and the total number of events 
seen. The number of observed events was consistent with the background estimates which 
dominate the SM WW/WZ signal. 

The trigger and offline electron identification efficiencies were estimated using Z — > ee 
events. The trigger efficiency was (98.1 ± 1.9)% |J. The electron identification efficiencies 
were found to be (74.5 ± 1.1)% in the central calorimeter and (61.9 ± 1.1)% in the end 
calorimeters. We studied the W — > jj efficiency for the jet cone size 1Z = 0.5 using the 
isajet and pythia |E6| event generators followed by detector simulation. The selection 
criteria were applied to these samples and it was found that the efficiency was ~ 50% for 
Pr < 250 GeV/c and that this decreased significantly for pj? > 250 GeV/c due to merging 
of the two jets into one. The efficiencies obtained from isajet were used to estimate the 
detection efficiencies of the WW {W Z) processes since they gave more conservative results. 

The overall event selection efficiency was calculated using the leading order event gen- 
erator of Ref. [17] to generate four-momenta for WW and WZ processes as a function of 
the coupling parameter values. A fast detector simulation was used to take into account 
the detector resolutions and efficiencies described above. Higher order QCD effects were 
approximated by a i^-factor = 1 + ^7ca s = 1.34 and a smearing of the transverse mo- 
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mentum of the diboson system according to the experimentally determined spectrum 
from the inclusive Z — > ee sample. The total selection efficiencies for the detection of 
SM WW and WZ events were estimated to be [14.7 ± 0.2 (stat) ± 1.2 (syst)]% and 
[14.6 ± 0.4 (stat) ± 1.1 (syst)]%, respectively. The systematic uncertainty (8%) includes: 
electron trigger and selection efficiencies (1%); $ T smearing and pr of the WW/WZ diboson 
system (5%); difference between the isajet and pythia programs for W — > jj efficiency 
parametrization (5%); statistical uncertainty for W — > jj efficiency parametrization (2%); 
and jet energy scale (3%). 

The expected signal for WW plus WZ production with SM couplings is 20.7±3.2 events 
based on the total integrated luminosity of 96.0 pb _1 . Figure [I] shows the distributions 
for candidate events from 1993-1995 data, total background estimate plus SM expectations, 
and SM expectations for WW and WZ production, after all selection criteria have been 
applied. There is no clear difference between the observed p^ spectrum and that expected 
from background plus SM WW and WZ prediction. 

Using the detection efficiencies for SM WW and WZ production and the background 
subtracted signal, and assuming the SM ratio of cross sections for WW and WZ production, 
we can set an upper limit at the 95% CL on the cross section a{pp — > W + W~X) of 76 pb. 

Since we observed no excess of high p^f events, large deviations from the SM trilinear 
coupling values are excluded. Limits on the anomalous coupling parameters were set by 
performing a binned likelihood fit to the observed pjf spectrum with the Monte Carlo signal 
prediction plus the estimated background. Unequal width bins were used to evenly distribute 
the observed events, especially those at the end of the spectrum. In each p^f bin for a given 
set of anomalous coupling parameters, we calculated the probability for the sum of the 
background estimate and Monte Carlo WW/WZ prediction to fluctuate to the observed 
number of events. The limits on the anomalous coupling parameters are from a combined 
likelihood fit to both data sets. The uncertainties in the background estimates, efficiencies, 
integrated luminosity, and higher order QCD corrections to the signal were convoluted with 
Gaussian distributions into the likelihood function. Uncertainties common to both analyses, 
e.g. theoretical uncertainties, were convoluted only once. 

In Fig. bounds on four pairs of coupling parameters are shown using A = 1.5 TeV. 
In each case all other couplings are fixed to their SM values. The one- and two-degree-of- 
freedom 95% CL contour limits (corresponding to likelihood function values 1.92 and 3.00 
units below the maximum, respectively) are shown as the inner curves, along with the S- 
matrix unitarity limits, shown as the outermost curves, which are obtained by evaluating all 
(i.e. WW, Wj, and WZ) processes. Figure 2(a) shows the contour limits when coupling 
parameters for WWj are assumed to be equal to those for WW Z. Figure 2(b) shows contour 



limits assuming HISZ relations [18]. In Fig. 2(c) and 2(d) SM WW^j couplings are assumed 
and the coupling limits for WWZ are shown. 

When SM WW'j couplings are assumed, the U(l) point {kz = 0, \z = 0, gf = 0) is 
excluded at the 99% CL. This is direct evidence for the existence of the WWZ couplings. 

Table II lists the 95% CL axis limits for three different values of A and four assumptions: 
(i) Ak = Ak 7 = Ak z , A = A 7 = \z', (ii) HISZ relations; (iii) SM WWj couplings; 
and (iv) SM WWZ couplings. The results with the SM WWj assumption are unique 
to WW/WZ production since the WWZ couplings are not accessible with Wj production. 
The results indicate that this analysis is more sensitive to WWZ couplings than to WWj 
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ones as expected from the larger overall couplings for WW Z than WWj 0. The dependence 
of the coupling parameters on A is clearly seen. Tighter limits are obtained when a larger 
value for A is used. When SM WW Z couplings are assumed, our limits on A 7 and Ak 7 with 
A = 2.0 TeV are not tight enough to lie within the S-matrix unitarity limit. 

In conclusion, we have presented limits on anomalous WW Z and WWj coupling param- 
eters which are the most stringent to date. They are significantly tighter than those from 
the analyses of the 1992-1993 data set [0,|J, and significantly better on An (but comparable 
on A) to those measured using W~f production with the 1992-1995 data set [[J. 

We thank U. Baur for useful discussions and D. Zeppenfeld for providing us with the 
WW and WZ Monte Carlo generators. We thank the staffs at Fermilab and collaborating 
institutions for their contributions to this work, and acknowledge support from the De- 
partment of Energy and National Science Foundation (U.S.A.), Commissariat a L'Energie 
Atomique (France), State Committee for Science and Technology and Ministry for Atomic 
Energy (Russia), CNPq (Brazil), Departments of Atomic Energy and Science and Education 
(India), Colciencias (Colombia), CONACyT (Mexico), Ministry of Education and KOSEF 
(Korea), CONICET and UBACyT (Argentina), and the A.P. Sloan Foundation. 
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TABLES 

TABLE I. Summary of signal and backgrounds. 



1992-1993 1993-1995 



Luminosity 


13.7 pb- 1 


82.3 pb" 1 


Backgrounds 








W+ > 2 jets 


62.2±13.0 


279.5±36.0 


QCD Multijet 


12.2L2.6 


104.3L12.3 


ti — > ev jjX 


0.9L0.1 


3.7L1.3 


Total Background 


75.3±13.3 


387.5±38.1 


Data 




84 


399 


SM WW+WZ prediction 


3.2±0.6 


17.5L3.0 


TABLE II. Axis limits at the 95% CL with various 


assumptions and three different A values. 


Couplings / A(TeV) 


1.0 


1.5 


2.0 


(i) A 7 = X z 


-0.42, 0.45 


-0.36, 0.39 


-0.33, 0.36 


Ak 7 = Akz 


-0.55, 0.79 


-0.47, 0.63 


-0.43, 0.59 


(ii) A 7 (HISZ) 


-0.42, 0.45 


-0.36, 0.38 


-0.34, 0.36 


Ak 7 (HISZ) 


-0.69, 1.04 


-0.56, 0.85 


-0.53, 0.78 


(iii) X z (SM WWj) 


-0.47, 0.51 


-0.40, 0.43 


-0.37, 0.40 


Ak z (SM WWi) 


-0.74, 0.99 


-0.60, 0.79 


-0.54, 0.72 


Agf (SM WWi) 


-0.75, 1.06 


-0.64, 0.89 


-0.60, 0.81 


(iv) A 7 (SM WWZ) 


-1.28, 1.33 


-1.21, 1.25 




Ak 7 (SM WWZ) 


-1.60, 2.03 


-1.38, 1.70 
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FIGURES 




50 100 150 200 2 

p e T (GeV/c) 



FIG. 1. pt distributions of the ev system for the 1993-1995 data set. The points with error 
bars represent the data. The solid histogram is the total background estimate plus the SM Monte 
Carlo predictions of WW and WZ production (shown as shaded histogram). The inset shows 
the predicted da/dp^ , folded with the detection efficiencies, for SM WW^f and WWZ couplings 
(lower curve), and for SM WW^j and the indicated anomalous WW Z couplings (upper curve). 
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FIG. 2. Limits on CP-conserving anomalous couplings parameters: (a) An 



Ak z 



\z', (b) HISZ relations; (c) and (d) SM WWj couplings. The inner and middle curves 



A = A- 

represent 95% CL one- and two-degree-of-freedom exclusion contours, respectively. The outermost 
curves show S-matrix unitarity bounds. A = 1.5 TeV is used for all four cases. The SM prediction 
is Ak = 0, A = 0, Agf = 0. 
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